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ABSTRACT 

The Rossi X-Ray Timing Explorer has made feasible for the first time the search for high-frequency (~100 Hz) 
periodic features in black hole candidate (BHC) systems. Such a feature, with a 67 Hz frequency, recently has 
been discovered in the BHC GRS 1915+105 (Morgan, Remillard, & Greiner). This feature is weak (rms 
variability ~0.3%-1.6%), stable in frequency (to within ~2 Hz) despite appreciable luminosity fluctuations, and 
narrow (quality factor Q ~ 20). Several of these properties are what one expects for a diskoseismic g-m°de m 
an accretion disk about a 10.6 A/© (nonrotating) to 36.3 M© (maximally rotating) black hole (if we are observing 
the fundamental-mode frequency). We explore this possibility by considering the expected luminosity modula- 
tion as well as possible excitation and growth mechanisms— including turbulent excitation, damping, a 
“negative” radiation damping. We conclude that a diskoseismic interpretation of the observations is viable. 

Subject headings: black hole physics — X-rays: stars 


1. INTRODUCTION 

In a series of previous papers (Nowak & Wagoner 1992, 
1993; Perez et al. 1997; hereafter NW92, NW93, and P97, 
respectively), we have discussed a class of modes in thin, 
Keplerian accretion disks that exist only in the presence of 
strong-field gravity (not in Newtonian gravity), some of whose 
eigenfrequencies depend primarily upon the mass and angular 
momentum of the black hole. The modes are trapped by 
general relativistic modification of the radial (k) and vertical 
(fU) epicyclic frequencies of free-particle, circular orbit per- 
turbations (P97, and references therein). 

These modes of oscillation are perturbations of the disk that 
are proportional to exp [/(oT + mcf))]. With disk angular 
velocity O(r), their corotating frequency is w(r) = or + mCl. 
Three classes of modes (/?-, g-, and c-modes) have been 
identified (P97, and references therein). However, here we 
shall be concerned solely with radial (m = 0) g-modes, which 
are the most relevant observationally. 4 These g-modes are 
trapped where to 2 < k 2 , in the region where k achieves its 
maximum value (at r — SGM/c~ for a nonrotating black hole) 
because, unlike in Newtonian gravity, k(t) rolls over at small r 
and vanishes at the inner disk edge, thereby creating a 
resonant cavity. The modes with the fewest radial nodes have 
relatively large radial extents, A r « GM/ c~. 

The frequencies (/ = — <t/ 27 t) of the radial g-modes are 

f = 714 F(a) (1 - «*) Hz, 
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(see P97). Here F(a ) is a known function of the dimensionless 
black hole angular momentum parameter a = cJ/GM 2 , rang- 
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4 These modes cover the largest area of the disk in a region in which the disk 
temperature is maximum (away from the uncertain conditions at the disk edge). 


ing from F( 0) = 1 to F(0.998) = 3.443. The properties of the 
disk enter only through the small correction term involving the 
disk thickness 2k, and the radial ( n ) and vertical (j) mode 
numbers, with 5 ~ 1 (as derived from the WKB solutions of 
P97). For thin-disk models in which h/r ~ 0.1L/L Ed d (in the 
mode-trapping region), e nj is typically on the order of 1% for 
the lowest modes (n ~ j ~ 1). Therefore, the mode frequency 
is relatively independent of disk luminosity. (The mode width, 
however, scales as Ar x c s 1/2 x h l/2 , where c s is the sound speed 
in the disk.) The radial g-mode is thus a good candidate for 
explaining the recently discovered high-frequency features in 
the black hole candidate GRS 1915+105. 

The observed frequency of this feature is 67 Hz, with a full 
width at half-maximum of 3-4 Hz (Morgan, Remillard, & 
Greiner 1996, 1997). This implies an effective Q = //A/ ~ 20. 
Despite factors of ~2 luminosity variations in the source, the 
frequency remained constant to within 1-2 Hz. The root mean 
square (rms) variability of the 67 Hz feature varied from 0.3% 
to 1.6% of the total observed X-ray luminosity, with the lower 
limit essentially being the Poisson noise limit. At its strongest, 
the rms variability is ~1% below 6 keV, and the rms variability 
is ~6% in the 12-25 keV band, which indicates that the 
spectrum of the feature is harder than the integrated source 
spectrum. (The spectrum of the source is fitted by a power law 
exponentially cut off at ~5 keV.) 

Below we discuss to what extent g-modes can explain these 
observations. In § 2, we further review the properties of these 
modes and make simple estimates of their luminosity modu- 
lation. In § 3, we discuss application to GRS 1915 + 105. In § 4, 
we discuss the role of turbulence and radiation as possiJ.e 
damping, growth, or excitation mechanisms. In § 5, we present 
our conclusions. 


2. LUMINOSITY MODULATION OF MODES 

For both the psuedo-Newtonian (NW92, NW93) and fully 
relativistic (P97) diskoseismology calculations, the equations 
are formulated in terms of the potential V — 6 / p, w 
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is the Eulerian variation of the pressure and p is the unper- 
turbed density. For modes with an azimuthal and time depen- 
dence ocexp [i(m<f> + or)], the potential V(r, 17, 4>, t ) can be 
separated into a radial, V r (r), and a vertical, V n (r}), component, 
where 77 = z/h{r) is the dimensionless vertical coordinate. The 
(small) fluid displacement vector, can be related to the 
Eulerian potential V(r, rj, <p, t) via the equations 

Ufa 4>, z, t) « (cu 2 - K 2 )' 1 

X exp + at)]VAdV r / dr), 

( 2 ) 

tz(r, <f>, z , t) « (a r - N 2 )~' 

X exp + at)]V r (dVJdz), 

where N : (r)) is the vertical buoyancy frequency, which we shall 
set to zero (see P97 for a discussion of the role of nonvanishing 
buoyancy frequency). We shall usually use A to denote a 
Lagrangian variation of a given quantity, where we follow the 
definitions of NW92, NW93, and P97. (For m = 0, g-modes 
have or = a 2 S k 2 .) 

The radial component of the fluid displacement approxi- 
mately satisfies the WKB relation 

d 2 W 

(o 2 c 2 (r, 0 ) — + (or - k 2 )(ot - Yfl 2 )JT = 0, (3) 

dr 

where Y(r) is a slowly varying separation function (akin to a 
separation constant in the WKB limit), c 5 (r, 0) is the speed of 
sound at the disk midplane, W = (k 2 - w 2 )~ ] dV r /dr, and = 
fl for a nonrotating black hole. We have taken the perturba- 
tions to be adiabatic, which means that the Lagrangian per- 
turbation of the pressure and density are given by A P/P = 
yA p/p « — y V • £. Utilizing the radial WKB approximation 
employed in equation (3), the (approximate) components of 
the divergence become 


« [(Yfl 2 - (i) 2 )/h 2 fl 2 io 2 ] exp (iat)V r (r)V v (ri), 


— « (ho>) exp (icrt)V r (r) . (4) 

dz dr) 

(We have also taken the unperturbed vertical barotropic 
structure to be in hydrostatic equilibrium.) 

For most simple a-models, the energy generation rate per 
unit volume is approximately aP(r, z)fl(r). The modes not only 
perturb the pressure (we ignore possible perturbations to a, 
and one can show that perturbations to (l are negligible), but 
they also perturb the locations of the disk boundaries. The 
variation of the luminosity is therefore 


f ro+triro) 
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r dr afLP(r, z) dz 


( 5 ) 


where P'(r\ z ') = P(r, z ) + AP(r, z), r ' = r + £(r), z' = 
z + £(z), and r h r Q , and Zb are the disk boundaries ( r 0 can be 


taken to go to oo without loss of generality). Transforming 
variables, the change in luminosity can be written as 

dr afl(/')r f dz P{r, z 
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where we have employed the above expression for the La- 
grangian pressure change and have employed the WKB ap- 
proximation in each harmonic term above. 

The luminosity fluctuation in general will have contributions 
from the fundamental frequency, as well as from the first two 
harmonics. We have calculated the luminosity modulation for 
modes in a disk with radial sound speed profile equal to that of 
a radiation pressure-dominated Shakura-Sunyaev disk with 
L/L Edd = 0.3. The exact radial profile does not greatly affect 
the estimates, as we can replace the quantity / P{r, O)H(r) dz 
with F(r), the disk energy flux, which is known from energy 
conservation. As the mode width A r ozh 1 ' 2 , the magnitude of h 
does affect the luminosity estimates. If we define the maximum 
magnitude (achieved at a radius r m ) of the vertical displace- 
ment vector £.(/•„„ z, t) to be =Wi(r m ), then for the case of the 
mode with one radial node in V, and two vertical nodes in 
V v (d&/dr) even about the midplane), 5 we have 


SL 

77 


(0.6%)^ sin (at + 0) - (9.8%)W 2 


X sin 2 (err + 0) + (Q.2%)^ 3 sin 3 (err -I- 0). (7) 


(Here we have used the psuedo-Newtonian approximation of 
NW92 and NW93 for our calculations.) This is the modulation 
of the bolometric luminosity. As the mode exists in the inner, 
hotter regions of the accretion disk, the modulation in re- 
stricted high-energy bandpasses will be greater. However, we 
have not properly accounted for radiative transfer effects. 
Dispersion in photon diffusion times over the extent of the 
mode will decrease the modulated luminosity for large optical 
depths. 


3. APPLICATION TO GRS 1915 + 105 

Although the above luminosity modulation is not large, 
Hi ~ 0.6-1. 6 produces, at the fundamental frequency, rms 
variability ~0.3%-1.6% (=[(0.6%)9C + (0.2%)^}/ \2), as 
is seen for the feature in GRS 1915+105. The first harmonic 
can produce this observed rms variability for M ~ 0.3- 0.7. 
(Therms variability is the coefficient of the sin 2 term divided 
by V8.) The modulation seen in GRS 1915+105 increases to 
6% if the energy band is restricted to the 12-25 keV range. 
This qualitatively agrees with what we expect for our modes. If 
we compare the bolometric luminosity variation of the mode 
to the bolometric luminosity from the region r « (6- 
20)GM/c 2 , the relative fractional variation increases to *6%. 

For the particular radial and vertical disk structure assumed 
here, the rms variability in the first harmonic is more signifi- 
cant than the rms variability in the fundamental frequency. As 
so little is known, both observationally and theoretically, about 

5 The first “even” mode about the midplane is not trapped in the region of 
the epicyclic frequency maximum. The first “odd” mode leads to a vanishing 
vertical integral in the above estimate of luminosity modulation. 
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disk vertical profiles, we cannot definitively state whether this 
feature is generic. Based upon energy arguments, however, we 
can say that generically disk modes can only produce rms 
variability ~0( 1%) if %€ ~ 1. We regard the question of which 
frequency, the fundamental or first harmonic, to associate with 
the observed 67 Hz feature as an open issue. 

Equation (1) indicates that a 67 Hz (fundamental) modula- 
tion can be produced in a disk around a black hole with a mass 
of 10.6 Mq for a nonrotating hole and 36.3 M e for a maximally 
rotating hole. There are few “natural” frequencies to invoke in 
a black hole system. If one were to appeal to a process 
occurring at the marginally stable orbit (r = 6GM/c 2 for a 
nonrotating black hole), not only would the luminosity from 
this region be lower than our mode estimates (owing to the 
no-torque boundary condition used in disk models), but also 
the required hole mass would increase by a factor of at least 
16/(3 N 3) ~ 3 (if the 67 Hz feature represents the fundamen- 
tal mode frequency). Our nonrotating hole, fundamental 
frequency, estimate of the required mass would yield an 
observed X-ray luminosity (Morgan et al. 1997) S30% of the 
Eddington luminosity, consistent with our assumptions above. 

4. MODE EXCITATION AND SELECTION EFFECTS 

4.1. Turbulent Damping and Excitation 

It is possible to use a parameterized stress tensor to estimate 
the effects of turbulent viscosity on the modes (NW92, NW93). 
The canonical energy of a radial mode is E c ~ o- 2 p(£: + £ 2 )AF, 
where AV is the volume occupied by the mode. Isotropic 
turbulence produces a rate of change dEJdt = ~EJ t, with 
r ~ \aor(h 2 /k; + h 2 /X:)\~' and A r , \ z , respectively, being the 
radial and vertical mode wavelengths. The corresponding 
quality factor is given by 

' Qj = (|cr|r)-‘ - Lr + (h/r)n 2 ]a, (8) 

as A ; ~ h/j and A r ~ (hr) ]/2 /n, where j and n are of the order 
of the number of vertical and radial nodes in any particular 
eigenfunction. Thus, for a 1, we can have high-mode Q. 

The other contribution to the fractional width Af/f of the 
corresponding feature in a power spectrum comes from the 
number of modes that are significantly excited. If we assume 
that this will include those whose wavelengths are greater than 
the maximum eddy size L ~ a U2 h divided by the corresponding 
Mach number ~a 1/2 , equation (1) yields 

A/// S A«(/i/r) ~ JhFr. (9) 

The minimum effective frequency width will be composed of 
this span of modes, each broadened by 1/0 ~ a. 

The above estimates are for isotropic viscosity. If the turbu- 
lence does not efficiently couple to the vertical gradients of the 
modes, then the mode 0-value is increased by a factor 
~UK/h) 2 (NW93). Aside from damping modes, turbulence 
can also potentially excite modes. Velocity perturbations in the 
disk, Sr, are made up of a mode component, Sr.u, and a 
turbulent component, &v T . Viscous damping arises from terms 
of the form 8v MJ 8v TJ , while mode excitation arises from terms 
of the form 8v Ti 8v rj (NW93). It is possible to make simple 
estimates of the magnitude of the turbulent excitation and to 
balance this against the turbulent damping (NW93). The modes 
are excited to an amplitude of |£~~| ~ a(h/\j i2 h, and ~a(A r /h) l/2 h, 
for isotropic and anisotropic viscosity, respectively (NW93). If 
turbulence is playing the dominant role in damping and exciting 


the modes, then we have the following constraints. For isotropic 
turbulence, the 0-value is large only for a 1; however, this 
implies a correspondingly small amplitude. For anisotropic vis- 
cosity, not only can we tolerate a larger a, but we also achieve a 
larger mode amplitude, although achieving the observationally 
required nonlinear mode amplitude is difficult in either case. 


4.2. Negative Radiative Damping 

As a first approximation, we took the modes to be adiabatic. 
In reality, we expect there to be small entropy changes due to 
various effects, the most notable one being radiative losses. If 
we have a radiation pressure- dominated atmosphere, as is 
likely in high-luminosity disks, one properly should use 


P P 

A P = y — Ap + y — As, (10) 

p s 


where s is the specific entropy. We can estimate 6 the effect of 
this term for a radiation pressure-dominated atmosphere, 
where one has 


P Ds c YP c 

4 — * v • — « — n 2 

S Dt K es p K ti 


( 11 ) 


( c is the speed of light, k £S is the electron scattering opacity, 
and we have invoked vertical hydrostatic equilibrium). 

For perturbations of scalars, the Lagrangian perturbation 
operator A commutes with total time derivatives D/Dt (see 
Lynden-Bell & Ostriker 1967). Using this fact, and the fact 
that for our modes DAs/Dt = iooAs, we can show that 


PAs 

s 


1 - i 


4 k„ s P to 
cH 2 


1 - i 


4q)T„/2 


AP - P£ r 


AP - P£ r 


d In D 2 
d In fl 2 


( 12 ) 


where T es is the electron scattering optical depth of the 
atmosphere. The term proportional to d In Q 2 /dr in the above 
is negligible in the WKB approximation. The Lagrangian 
variation of the entropy is therefore directly proportional to 
the Lagrangian variation of the pressure. 

We do not know the true radial and vertical structure of a 
disk model that correctly describes the spectral observations of 
GRS 1915 + 105. However, for a “standard” Shakura-Sunyaev 
a-disk model, one can show, with to ~ k, that cot es h/c ~ a~ l 
and hence is likely to be ^>1. Taking this to be the case, we can 
combine equations (10) and (12) to yield AP/P = yAp/p, 
where 

y’ * 7^1 + i h j = 7(1 + «*')• (13) 

We have subsumed our ignorance of the disk’s vertical struc- 
ture into the parameter a', which we expect to be of O(a). 

To estimate the effects that the nonadiabatic terms have 
upon our modes, we can substitute the above into equation (3) 
with d 2 W/dr 2 « -k 2 W and then expand it about the adiabatic 
solution. Specifically, we then have a dispersion relation that 
can be written in terms of a function G(co) = 0, to which we 
are adding a function H(co) due to radiation damping. If c^ is 

6 Throughout this section, our estimates will be based upon vertically 
averaged quantities rather than fully 2 -dependent quantities. 
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the unperturbed mode frequency and 6o> is its perturbation, 
we then have to first order 


dG 

do) 


So + H(c j () ) = 0. 


The mode O-value then becomes 
. dG/d In wL 


Q 1 So 1 


2Y / a \ 2 


H(co 0 ) 


a' \cX) ’ 


(14) 


(15) 


where, on the right-hand side of the above, we have used 
equation (3) and have taken the term YH 1 to be dominant 
(i.e., 1 5 1). The minus sign indicates mode growth', that is, 
the modes are unstable to radiative losses. As typically k~ x > h, 
Y > 1, and or < 1, it is fairly easy to obtain \Q\ :» 1. 


4.3. Selection Effects 

As discussed in § 2, the modes must approach nonlinearity 
in order to be observationally relevant. If isotropic turbulent 
excitation determines the mode amplitude, then there is a 
natural maximum vertical perturbation amplitude with 
1. If the turbulence is anisotropic, on the other hand, then it is 
possible to achieve W - 1 , so long as a ~ (h/\ r ) 1/2 £ 1 . For 
this case, we expect Q-values ~af 3 . Radiative excitation also 
naturally leads to ~ 1. At some point, currently unknown, 
damping effects and nonlinearities must limit this radiative 
mode growth. As a zeroth-order approximation, we take the 
linear growth rates calculated above to be relevant to the 
near-nonlinear regime and thus to be equal to the saturated 
damping rate. The above radiative (7-value then gives a rough 
estimate of the width of observationally relevant modes. 

The low rms variability of the modes means that only modes 
with high (7-values w ill be detectable. Observational limits for 
Q ~ 20 were rms 0.3% and detectable rms oc Q~ l/ 2 for 
narrow modes. Considering the dominant contributions as n 
and; increase, turbulence leads to Q oc j~ 2 (isotropic) and Q oc 
n ^(anisotropic), whereas radiative damping leads to Q oc 
j /n . The limiting detectable rms variability is therefore linear 
in n and;. Turbulent excitation leads only to high (7-values for 
modes with low; (isotropic) or low n (anisotropic). Radiative 
damping, on the other hand, leads only to high (7-values for 
modes with low n and high;. We note, however, that the mode 
calculations require a smooth, well-behaved background to 
perturb. If small-scale turbulence truly is responsible for disk 
viscosity, our unperturbed hydrodynamic equations are rele- 
vant only for the largest scales, again favoring both low n and 
low;. For all cases considered above, a high (7-va lue is most 
easily achieved for a « 1. The 67 Hz mode seen in GRS 
1915+105, if a diskoseismic mode, therefore is likely one with 
low n and ;, and a <&; 1. 


5. CONCLUSIONS 

The main motivations for attributing the 67 Hz feature in 
GRS 1915 + 105 to a diskoseismic mode are that (1) these 
modes are related to a “natural” frequency in the disk (i.e., the 
maximum epicyclic frequency); (2) their spectra are expected 
to be characteristic of the inner, hottest regions of the disk; (3) 
their frequencies are relatively insensitive to changes in lumi- 
nosity; and (4) they have low rms variability. This latter 
feature, although in agreement with the observations, is the 
strongest constraint. These modes cannot be applied to sys- 
tems that show £10% rms variability over a wide range of 
energy bands. We have identified two mechanisms, turbulent 
excitation and negative radiation damping, that naturally lead 
to appreciable mode amplitudes with high (7-value for a <&: 1. 

Again, we stress that the only other “natural” disk fre- 
quency, the Keplerian rotation frequency at the last marginally 
stable circular orbit, would require a black hole mass a factor 
£3 larger. Furthermore, less flux is emitted from that region 
than from the g-mode region. 

We saw that there were two requirements for the g-modes to 
be observationally detectable. First, the g-modes had to be 
wide, or equivalently the disk had to be luminous, assuming 
that cf 2 oc V 12 . Second, the g-mode amplitudes had to 
approach the nonlinear regime. Ideally, one should perform 
numerical simulations of the nonlinear development of the 
g-modes to determine mode widths, (7, and luminosity modu- 
lation more accurately. It is interesting to note that the MHD 
simulations of Stone et al. (1996) do show copious p-mode 
production associated with the turbulence (Gammie 1996). 
(MHD simulations have yet to be performed for rotation and 
epicyclic frequency profiles relevant to the g-mode-trapping 
region.) 

These g-mode oscillations should also exist in accretion 
disks around compact (soft equation of state), weakly magne- 
tized ( B < 10 8 G) neutron stars. Under these conditions, the 
inner radius of the disk will be less than those radii where the 
g-mode exists. Again, large-amplitude (rms » 1%) features 
cannot be explained with these modes. 

Even if the 67 Hz feature seen in GRS 1915+105 does not 
turn out to be a diskoseismic mode, it points out two impor- 
tant lessons. First, black hole candidate systems can pro- 
duce relatively stable, high-frequency features. Second, the 
Rossi X-Ray Timing Explorer is capable of detecting and 
characterizing these features despite their weak variability. 
The search for diskoseismic modes in this system and other 
BHC has therefore become a viable and worthwhile pursuit. 
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